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ABSTRACT 

We detect the HD 32297 debris disk in scattered light at 1.6 and 2.05 /.im. We use these new observations 
together with a previous scattered light image of the disk at 1 . 1 to examine the structure and scattering 
efficiency of the disk as a function of wavelength. In addition to surface brightness asymmetries and a warped 
morphology beyond ~1"5 for one lobe of the disk, we find that there exists an asymmetry in the spectral 
features of the grains between the northeastern and southwestern lobes. The mostly neutral color of the disk 
lobes imply roughly 1 /im-sized grains are responsible for the scattering. We find that the asymmetries in color 
and morphology can plausibly be explained by HD 32297's motion into a dense ISM cloud at a relative velocity 
of 15 km s -1 . We model the interaction of dust grains with HI gas in the cloud. We argue that supersonic ballistic 
drag can explain the morphology of the debris disks of HD 32297, HD 15115, and HD 61005. 
Subject headings: stars:individual(HD 32297) — stars :individual(HD 15115) — stars :individual(HD 61005) — 
circumstellar disks — methods:n-body simulations — ISM 
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1. INTRODUCTION 

Debris disks are created by the outgassing and collisions of 
planetesimals that may resemble analogues to the comets, as- 
teroids, and dwarf planets in our Solar System. Debris disks 
show us the conditions and environments of planet building 
for stars different than the Sun through their structure and 
mineralogy. 

As more disks are resolved at multiple wavelengths, they 
reveal their intrinsic complexity. Both the morphology and 
spectra of disks change in ways that are not fully understood. 
In order to probe the intrinsic properties of a disk, effects that 
are due to the surrounding environment of a star must be ac- 
counted for. 

HD 32297, an A5 star at 112±11 pc with a recently dis- 
covered edge-on debris disk, appears to be a promising ob- 
j ect of study. Fi r st noted as a star wi th an IRAS excess 
(ISilverstonel2000l) . ISchneider et al.l(l2005l Hereafter S05) dis- 
covered an edge-on disk structure in scattered light at 1 . 1 /im. 
The disk extended out to ^400 AU and showed surface 
brightness asymmetries between the two lobes of the disk. 
This morphology appears to be common among known edge- 
on disks like /3-Pict oris and HD 15115, which both show 
similar asymmetries dKalas & Jewit t 1995; Golim owski et al.l 
I2006t iKalas et al.l l2007t iDebes et al.l l2008alh The total 
L/ R /L*=0.27% makes the disk optically thin S05. 

Observations of HD 32297's outer disk in the R band re- 
vealed that it may be interacting with the surrounding ISM. 
The disk app eared signifi cantly warped at larger separations 
from the star dKalasll2005l hereafter K05). The NE half of the 
disk appeared to bow away from the midplane of the disk. An 
extrapolated surface brightness profile implied blue colors for 
the disk which K05 interpreted as resulting from sub-micron 
silicate dust grains. K05 suggested that the bowed structure 
results from an ISM-disk interaction. 

The disk was also spatially resolved at 12 and 18Atm 
in ground based observations with Gemini ( Mo erchen et al.l 
120071: iFitzgerald et a l. 2007b). There is a suggestion that at 
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distances >0."75 the disk is brighter on the NE side, i.e. re- 
yersed from the asym metry seen in scattered light images 
dMoerchen et al. 2007). Models of the thermal emission and 
brightness temperature curves as a function of distance based 
on the 12/im/18/!m colors suggest that the disk has a ring- 
like structure with an interior clearing out to ~60-80 AU 
( IFitzgerald et alJl2007bT:lMoerchen et ai]|2007l) . 

Finally, HD 32297's disk has been resolved with CARMA 
at 1.3 mm. The CARMA image shows a marginally resolved 
asymmetry in the disk that appears to match asymmetries seen 
in scattered light (Maness et al. 2008). Furthermore, the emis- 
sion at 1.3 mm is well in excess of that inferred by shorter 
wavelength thermal emission. Models of the disk so far rely 
on analytic models of dust grain optical behavior representa- 
tive of "dirty ice" mixtures, but there are few constraints on 
the true composition of the dust gra ins in the disk; the pre- 
sumption of water ice is premature (Fitzge rald et al J [2007b; 
iManess et alj|2008l) . 

In addition to dust, a small amount of gas is present in the 
HD 32297 system. Based on the detection of narro w circum- 
stellar Na I absorption in the spectrum of HD 32297. iRedfieldl 
(120071) determined that at most 0.3M ffi of gas resides within 
the disk. The CARM A observations also provided an upper 
limit of 95 M e of gas dManess et alj|2008l) . 

In this article we present two new NICMOS coronagraphic 
scattered light images of the disk at 1.6 and 2.05/im. We ex- 
amine how the morphology of the disk changes with wave- 
length and how that might lead to new constraints on the 
composition of the dust in the disk. The properties of the 
disk suggest that it is interacting with the surrounding ISM. 
We model this interaction and explore how it affects the ob- 
served properties of the disk. In ^2] we present the observa- 
tions of HD 32297, in |3] we measure the morphology, surface 
brightness, and scattering efficiency of the dust as a function 
of wavelength. Finally, we present models of the interaction 
between the ISM and HD 32297, HD 15115, and HD 61005 
ing] 

2. OBSERVATIONS 

We observed HD 32297 with the Hubble Space Telescope's 
(HST) Near-Infrared and Mu lti-Object Spectrometer (NIC- 
MOS) dThompson et all 1998) using the coronagraph in Cam- 
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era 2. We used the F160W and F205W filters with cen- 
tral wavelengths of 1.6 and 2.05^tm. Additionally, we re- 
analyzed archival data of HD 32297 in the Fl 10W filter taken 
by S05 All new observations were taken on 3 October 2007 
for HD 32297. 

In between visits to HD 32297 we observed HD 31489, the 
point-spread function (PSF) reference. Accurate disk pho- 
tometry requires a well constrained scaling ratio between the 
target star and its PSF reference. To achieve the most accurate 
scaling, over the three visits we obtained three short 0.6s ex- 
posures of both the target and PSF reference star without the 
coronagraph in both filters. 

In F160W, the reference saturated in each 0.6s exposure. 
To determine the proper scaling factor for F160W, we used an 
annulus centered on each star with a radius between 6 to 21 
pixels to determine the total counts/s measured. In F205W, 
this was not an issue for either star, and so we used a simple 
circular aperture with a width of 21 pixels. In F160W we 
found the ratio of flux densities for HD 32297 to HD 31489 
to be 0.54±0.02. In F205W, we found a ratio of 0.53±0.01. 

We also used the archival observations of HD 32297 and 
HD 83870. HD 83870 was reported in S05 to yield the best 
PSF subtraction using pixel offsets of [+0.304,+0.143] and a 
ratio of HD 32297 to HD 83870 of 0.440±0.005. Short ex- 
posures of both stars were also available, though both were 
saturated at their cores. We re-determined the scaling using 
the same procedure as we used for the F160W filter and ob- 
tained a ratio of 0.441 ±0.005, in excellent agreement with 
that determined by S05. 

In order to determine the best subtraction we minimized 
a chi-squared metric on a region of the target image dom- 
inated by the star's diffraction spikes. We assumed that 
good subtraction of the diffraction spikes corresponded to 
the best subtraction o f the PSF within the region of interest 
dSchneider et al.l 120011) . We iteratively created subtractions 
for combinations of scaling and pixel offsets until we found 
an image that produced the lowest chi-squared measure. We 
searched within l-a of the scaling ratios we determined above 
and within ±1 pixel to find the best x and y pixel offsets. 

To quantify the systematic effects on the photometry, we 
repeated the subtractions varying the PSF scalings and offsets 
by ± 1 a from the minimum chi-square solution found above. 
Using a rectangular photometric aperture matched to the size 
of the disk, we found the standard deviation in the disk flux 
densities from this suite of subtractions. We then propagated 
this uncertainty into the total uncertainty in the flux density of 
the disk per pixel. 

We observed HD 32297 at two telescope orientations dif- 
fering by 24°. This is essentially an azimuthal dither that al- 
lows astronomical objects that rotate with the field of view 
to be distinguished from rotationally invariant instrumental 
artifacts. The two differentially-oriented PSF-subtracted im- 
ages in each filter band were geometrically rectified and cor- 
rected for the linear optical distortion (X JC£I / f ,= 75.95 mas 
pixel -1 , Y sca / f =75.42 mas pixel -1 ) at the Camera 2 corona- 
graphic focus. The geometrically corrected image pairs were 
rotated about the location of the occulted target (as determined 
through the target acquisition process) and averaged together. 
The final, combined images are shown in FigureQ] 

In practice, instrument artifacts are rarely constant in time, 
so subtraction residuals remain and can contaminate photo- 
metric measurements of a disk. We therefore checked that the 
total disk flux measurement at each telescope orientation was 
consistent to the above quoted uncertainties. 



3. RESULTS 

3.1. Disk Geometry 

To measure any warping of the HD32297 disk system, we 
calculated the disk position angle (PA) as a function of radial 
distance. To calculate the PA, we performed the following 
analysis. First, we separated the images at each spacecraft 
orientation for all three filters to create six independent im- 
ages of the disk. In each image, we highlighted the midplane 
of the disk by performing an unsharp mask; we subtracted a 
Gaussian-smoothed disk image from the original image. Then 
we used a rectangular aperture of width 3 pixels in the radial 
direction and height 8 pixels in the direction perpendicular to 
the midplane of the disk over the disk image. We centered 
the aperture on the disk midplane at radial distances of 6 to 
40 pixels. In each aperture we summed the surface bright- 
ness along the 3 pixels in the radial direction radial direction 
and found the center of the disk by fitting the resulting cross- 
section with a Gaussian. For each radial distance we took a 
median of the six independent measures of the PA from each 
image and calculated an uncertainty based on the standard de- 
viation of the measures. The NE lobe has a median PA of 
46.3±1.3 (i.e., 226° ± 1° and the SW lobe has a median PA 
of 229° ± 2°. Our PA measurement differs from that reported 
in S05. However, the Fl 10W PA was corrected to 47.6° ± 1 
(i.e., 227.6° ± 1) in K05 by private communication from G. 
Schneider. 

The PA as a function of radial distance as measured from 
the stellar center is shown in Figure [2] The S W lobe of the 
disk shows a significant non-zero slope to the PA as a function 
of distance from the star, i.e., the disk midplane is curved. The 
NE lobe shows that the PA is uniform to within the uncertain- 
ties. However, visual inspection of the outer edge isophotes 
of the NE lobe with respect to the average PA measurements 
suggests that there might be some slight warping. 

To determine the extent to the warping of the SW lobe, we 
performed a least squares linear fit to the PA as a function of 
radial distance from the star. We found a slope of 2.1 ±0.4°/". 
K05 found significantly different PAs for both lobes compared 
to S05; they reported a PA of 245° ± 2° for the SW lobe but 
they did not specify at what distance this was measured. How- 
ever, if we extrapolate the curvature we see out to 8", the ra- 
dial midpoint of where the disk was detected in R, we would 
expect the PA to be 246.1°±3.4°, consistent with what was 
reported in K05. 

3.2. Surface Brightness Profiles 

To study the surface brightness profiles of the HD 32297 
disk we placed 3x3 pixel apertures on the images centered 
on the disk midplane. We measured the surface brightness 
of the disk in each aperture as a function of distance from 
the star and wavelength for both disk lobes. We used back- 
ground apertures placed just beyond the disk to estimate the 
level of any systematic uncertainties. These background aper- 
tures were located at similar distances from the star as the 
photometric apertures. 

Figure [3] shows the profiles in each filter. Overall the pro- 
files have a qualitatively similar shape. Of particular interest 
is the amount of asymmetry between the lobe surface bright- 
nesses; S05 noted that the NE lobe of the disk was dimmer 
than the SW lobe beyond l."7. This asymmetry is confirmed 
at new wavelengths by our data. We find that the asymme- 
try is more prominent at F160W and F205W, extending from 
about l"to2"5. 



HD 32297 at multiple wavelengths 



3 



The surface brightness profiles in both lobes and at all 
wavelengths cannot be described by a single power-law, rather 
they are best fit with broken power-laws. We tabulate the best 
fitting parameters for the power-laws in TableQ]along with the 
radius of the break (Rbreak), the power-law index before Rbreak 
a and the power-law index exterior to Rbreak- The radius at 
which the break occurs is slightly different between the two 
lobes, occuring closer to 90 AU on the NE lobe and 110 AU 
on the SE lobe of the disk, though this is at the limit of our 
spatial resolution and may not be significant. Both lobes are 
characterized by a slope less than r" 2 in the interior of the disk 
indicating a rise in the optical depth of scatterers, followed 
by a drop-off in surface brightness of ~r~ 3 beyond the break 
corresponding to a slow drop-off in optical depth beyond the 
break. 

3.3. Scattering Efficiency vs. Wavelength 

The next step in our analysis is to infer some scattering 
properties for the dust grains in orbit around HD 32297. 
While resolved thermal images of the disk exist, it is not clear 
if any of the observations to date constrain the composition 
of the material in the d isk (Silverstone 2000; Moerche n et al.l 
I2007t iFitzgerald etall l2007bt iManess etal.1 120081) . The 
SED at 12/xm shows no indication of emission around the 
strong 10/im silicate line, typical o f colder debris disks , 
or debris disks with la rger particles (Moe rchen et al.l 120071: 
IFitzgerald etaill2007bl) . In HD 32297's case, the tempera- 
ture of the dust is estimated to be ~170-180 K, hot enough 
to produce a silicate line if small silicate grain s are present 
dMoerchen et al.l 120071: IFitzgerald etail l2007bl) . For these 
types of debris disks in particular, scattered light images pro- 
vide one of the few avenues with which to infer dust compo- 
sition. 

Translating observed surface brightness profiles into infor- 
mation about the dust grains present requires removing the 
spectral contribution from the star itself to reveal the spectral 
character of the disk. In practical terms, the scattering effi- 
ciency of dust grains, Q sca , is proportional to the ratio of the 
surface brightness of the disk (SB„) to the flux density of the 
central star (F„). We will measure Q sca by dividing our im- 
ages (SB„) by the F^. 

For HD 32297, the direct imaging in each filter allows us 
to accurately estimate the flux density of HD 32297 for the 
different wavelengths observed. Since HD 32297 was not sat- 
urated in the F160W and F205W filters for our 0.6s exposures, 
we found that F„=0.93 ±0.02 Jy and 0.61 ±0.01 Jy forF160W 
and F205W respectively. 

The Fl 10W images pose a slight problem for estimating Fj, 
since HD 32297 saturated the detector for the 0.6s exposures. 
However, we can use the wings of the PSF to infer the total 
flux density of the star. We modeled the HD 32297 PSF with 
the tinytim software package (V6.2) 4 , assuming an A5 star. 
We used the same annulus described in § for determining our 
PSF scalings to calculate an aperture correction in our syn- 
thetic images between the flux in the wings to the total flux in 
the PSF. We assumed a 5% error in the multiplicative correc- 
tion term, or 9.6±0.5. Our resulting flux density in F110W 
is 1.55±0.08 Jy. S05 estimated the F110W flux density of 
HD 32297 using SYNPHOT. They estimated a flux density of 
1 .46±0.04 Jy, consistent with our measurement uncertainties. 

We used our 3x3 pixel apertures as described in ^3.21 to 
convert our surface brightness measurements into scattering 

4 http://www.stsci.edu/software/tinytim/tinytim.html 



efficiencies for distances from the star ranging from 0"6 
(67 AU) to 3" (337 AU). Figures H and |5] show the scatter- 
ing efficiency as a function of wavelength for 76, 127, and 
313 AU, showing a representative range of colors for the disk 
as a function of distance. In particular, both lobes of the disk 
are roughly neutral in the interior of the disk and blue at the 
outer edges. The SW lobe has a red zone between 127 AU 
and 313 AU. 

The surface brightness asymmetries appear to be the result 
of a color difference between the lobes. The color change be- 
gins at ^127 AU, corresponding roughly to where the warp in 
the Southwestern lobe becomes noticeable. The mostly neu- 
tral colors in the near-IR favor grains in the disk that are close 
in size to (or larger than) the wavelength of light observed, or 
r^l/im. Both silicates and water ice would be primary candi- 
dates to produce the observed efficiencies. 

4. ASYMMETRIES AND WARPS DUE TO INTERACTION WITH THE 

ISM 

The discovery of <5-Velorum 's high velocity inte raction with 
the ISM in the Local Bubble (Gas paret all l2008) and the pe- 
culiar morphology of the debris disk HD 61005 demonstrate 
the need to also cons ider the interaction of ISM gas with cir- 
cumstellar material dHines et alJ I2007L hereafter H07). A 
general picture of how t he ISM may affect circumstella r disks 
is nicely summarized in lArtvmowicz & Clampinl d!997l here- 
after AC97) as well as in H07 and references therein. We 
argue that the simplest explanation for the color and bright- 
ness asymmetries, as well as the perturbed morphology of the 
outer disk, is the interaction of HD 32297 with a cloud of ISM 
material as first suggested by K05. 

The behavior of HD 32297's surface brightness profiles are 
similar in behavior to the disks of AU Microscopii and j3- 
Pictoris, both of which show breaks in their radial bright- 
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ness profiles dHeap et al.ll2000t iGo limowski et all 120061: 



Liu 



2004; Kalas et all [20041: iMetchev et alJ 120051: iGraham etal 



2007; Fitzg erald et alJ2007al) . These breaks are interpreted as 
arising from the "birth-ring" of large planetesimals that feed 
both the inner dust population through Poynting-Robertson 
or wind drag and the outer dus t population through rad iation 
pressure and mutual collisions ( Stru bbe & Ch iang 2006J). The 
outer disk is then composed of increasingly weakly bound 
dust grains, or grains that are unbound and leaving the sys- 
tem (so-called j3 meteoroids). 

There are two components to an ISM cloud that dust grains 
both bound and unbound in a disk will interact with. First is 
the hydrogen gas which makes up the bulk of the mass in a 
cloud, and thus is the dominant perturbation on the dynam- 
ics of the grains through supersonic ballistic gas drag. The 
second component is the small ISM grains in the cloud which 
will primarily affect the composition and grain size distribu- 
tion of the disk grains through collisions with dust in the disk. 
ISM grains are repelled by radiation pressure of the central 
star at some radius, so this effect is primarily in the outer re- 
gions of a debris disk (AC97). Both of these components can 
have an observable affect on the morphology of the disk. In 
this Section we develop for the first time a simple model for 
a disk interacting with the gaseous part of the ISM. We apply 
it to HD 32297 and two other disks that have peculiar mor- 
phologies, HD 15115 and HD 61005. 

4.1. 77ie Interaction of Dust Grains with ISM Gas 

We consider the situation where a collisional cascade has 
resulted in a steady flux of unbound (3 meteoroids traveling at 
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escape speed from the system and encountering the ISM gas. 
Since the relative velocity of the dust grains is on the order 
or greater than the sound speed of the gas up to a gas tem- 
perature of —8000 K (the maximum temperature for a neutral 
cloud), the interaction of a single dust grain with the HI gas 
in the neutral ISM is akin to supersonic gas drag felt by mete- 
orit es entering the Earth's atmosphere with a drag coefficient 
D g dPecina & Ceplechalll983t lHoppdll937l) . In this model, 
the acceleration (a) a single dust grain feels from the grav- 
ity of the central star, the radiation pressure from the central 
star, and the drag fro m the ISM gas at a distance r is (e.g. 
Nak agawa et al]l 19861) 



GM^l-P) 



r-Dg (v dus t-Vg as ) 2 



(1) 



where the drag c oefficient is dHoppd 119371; 
iPecina & Ceplechall 19831) 
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Here, A is a shape parameter and T is the drag coefficient for 
a supersonic spherical particle {AT — 1), p gas and pd us , are the 
densities of the ISM gas and the disk dust grains respectively, 
Vdust and Vgas are the velocities of a dust grain and the ISM 
flow respectively. We define (3 as the ratio of the force of 
radiation pressure to the gravitational force exerted by the star. 

An order of magnitude estimate for the radial distance at 
which the dust grains are significantly perturbed from their 
original trajectories can be obtained by setting a=0 and as- 
suming weakly bound grains, i.e. (3 —0.5, in EquationQ] 
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where v re i is the relative velocity between the disk and the 
ISM. The assumption of (3 ~ 0.5 implies a particular grain 
mass and density for a given dust grain composition. (3 for 
a perfectly abso rbing spherical dust grain with radius rdust is 
(Burn s et al.lfl979l) 
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We can investigate /^deflect as a function of /9 gas , v re i, and L* 
for a star moving through a typical ISM cloud. If we combine 
Equations [3] and [4] for (3=0.5 we remove the dependence on 
both dust grain size and density and create a scaling law for 

^deflect^ 
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To synthesize model debris disk images to compare to our 
data, we integrated Equation Q] for various dust and cloud pa- 
rameters using a modified N-body integrator. We released 
thousands of particles with initial velocities determined by the 
birth ring orbital parameters: with an initial semi-major axis, 
uniformly distributed inclinations from 0° to 5°, uniformly 
distributed eccentricities <0.1, uniformly distributed longi- 
tudes of ascending node, uniformly distributed longitudes of 
pericenter, and true longitudes uniformly distributed between 
and 2tt. We assumed the dust grains have /3=0.51, that is, 
they are about the blow-out size for the star. A more thorough 



model could self-consistently calculate f3 for a given compo- 
sition and grain size. This approach is unjustified for most de- 
bris disks as there are currently few compositional constraints. 

We calculated each particle's trajectory with a Bulirsh- 
Stoer integrator (Pre ss et aill 19921) with an initial timestep of 
5 years, for a total of 5000 yr. We recorded particle positions 
every 5 years and accumulated them in a two dimensional his- 
togram of projected particle density 

We constructed a final disk image by calculating the 
surface brightnesses of the columns of particles in each 
pixel. We then Gaussian smoothed the resulting image to 
model convolution with the PSF of a telescope. For ex- 
ample, a simulation of a NICMOS image of HD 32297 
at F160W has a pixel size of 75 mas (corresponding to 
8.5 AU at HD 32297's distance) and the image uses a Gaus- 
sian smoothing kernel with a FWHM of ~137 mas. This 
method of using histograms to represent the steady-state dis- 
tribution of_duj>t_cj£udshave_be^ other 
authors JWilner et al.1 [20021: iMoro-Martfn & Malhotrall200"l 
iMoran et al.ll2004HStark & Kuchnerl2008h . 

4. 1. 1. An ISM Gas Interaction Model for HD 32297 

Modeling the structure of the HD 32297 disk ab initio 
would require knowing the mass, r e //,and luminosity of the 
star, the relative motion of the cloud compared to the disk, 
the composition of the disk grains, and the density of the gas 
cloud, all of which are poorly constrained or unknown. There- 
fore, we have simply chosen initial conditions to qualitatively 
match the observed morphology. With this approach we can 
1) demonstrate that this scenario can plausibly explain many 
features within the disk using reasonable initial conditions, 
and 2) back out some weak constraints on the disk grains and 
cloud parameters that appear to fit the data. 

Before we experimented with the dynamical models, we 
re-examined the stellar parameters of HD 32297. Litera- 
ture values for the sp ectral ty pe to HD 32297 range from AO 
dPerrvmanet al.ll 19971) to A5 (Fitzgerald et"afll2007bl and ref- 
erences therein). We conducted a least-squares fit of Kurucz 
models to HD 32297's B VJHK, photometry to determine both 
T e ff and luminosity. The model that best fits the photometry 
has T e fj =7750 K and luminosity = 5.7 L©. 

However, thes e values place HD 3 229 7 below the ZAM S on 
the tracks of both lSiess et al.ld2000l) and lD'Anton"aldl987l) . i.e. 
in an unphysical regime. The 1-er uncertainty o n the distance 
to the star in the new Hipparcos reduction of Ivan Leeuwenl 
(2007) is 12 pc. Keeping the star at 7750 K and moving it 
farther away by la, i.e. from 1 12 pc to 133 pc, would increase 
the inferred luminosity to 8. 1 L©; this is the luminosity on the 
ZAMS for a star of 1.7 M© on the Siess tracks. We take these 
values to be the most likely stellar luminosity and mass. 

One caveat to the above analysis is the possibility of extinc- 
tion from the intervening ISM or circumstellar material. If the 
intrinsic T e ff of the star is closer to —10000 K, as would be 
appropriate for its earliest literature spectral type of AO, then 
considerable extinction would be necessary to fit the photom- 
etry. An ISM-like A v =0.72 also provides a good fit to the ob- 
served BVJHKs photometry for T e// =10000K and log g=4.5 
Kurucz models. In this case the luminosity at the nominal 
distance is 13.8 L©, which is also well below the ZAMS at 
that T eff . At T e// =10000K, the ZAMS has luminosity -30 
L© for a M— 2.3 M©. That would require moving the star 
out to 165 AU, or 4.4<r on its Hipparcos dist ance. Given th e 
low ISM Na column de nsity toward the star (Redfield 2007), 
Fitzge rald et alj d2007bl) dismissed the possibility of high A v , 
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but other possible A v , T e ff combinations could result in an 
allowed luminosity. 

Alternatively, the star could experience gray extinction 
from circumstellar material to explain the discrepancy be- 
tween the inferred luminosity and the ZAMS luminosity. 
This, however, would imply an optically thick circumstel- 
lar disk comprised of fairly large dust grains. This situation 
seems unlikely given the low L«/L* measured which is mor e 
in line with optically thin dust disks (Fitzger ald et al.| [2007b). 
Also, the inferred amount of circumstellar gas is small, argu- 
ing for a late stage of evolution for t he disk where r elatively 
little primordial gas and dust remain (Redfield 2007). 

In addition to the stellar parameters, we must consider the 
relative velocity vector between the disk and the ISM flow, 
j3, and D g to fully determine the initial conditions for our 
simulations. The relative velocity vector can be weakly con- 
strained by the local dispersion of ISM cloud velocities and 
HD 32297's proper motion. The other parameters are more 
difficult to constrain. 

In order to determine a v re i consistent with both the ob- 
served proper motion and radial velocity of HD 32297 as well 
as the morphology of the warped disk, one must make some 
assumptions about the impacting cloud's velocity vector rela- 
tive to that of the star. From observations of the circumstellar 
and interstellar Na lines in HD 32297's spectrum the radial ve- 
locity of HD 32297 is r oughly 20 km/s while the ISM's radial 
velocity is ~ 24 km/s ( Red field! l2007h . HD 32297's proper 
motion is primarily to the South at a speed of 15 km/s (assum- 
ing a distance closer to 133 pc). If the cloud only moves in 
the radial direction, the magnitude of the relative velocity be- 
tween the cloud and HD 32297 would be 16 km/s. The cloud 
may have additional proper motion relative to HD 32297. 
If we assume that the cloud's velocity vector is parallel to 
HD 32297's motion on the sky that would produce a mini- 
mum relative velocity, while if it was moving anti-parallel it 
would produce the maximum relative velocity. Based on the 
dispersion in velocity of nearby ISM clouds, we assume the 
magnitude of th e cloud's velocity relative to that of HD 32297 
is ~14 km s" 1 (Red field & Linsky||2008l) . the minimum and 
maximum relative velocities between the cloud and HD 32297 
would be 4 km s" 1 and 29 km s" 1 . 

We varied the relative velocity vector and D g until we found 
qualitatively similar morphologies between our models and 
HD 32297. The best morphological match came from a rela- 
tive velocity vector with a magnitude of 15 km/s at an inclina- 
tion of 55° with respect to the y-z plane in coordinates where 
the x-axis is aligned with the disk midplane and the z axis is 
perpendicular to the midplane. This represents parameters in 
the middle of our estimates above. 

For this model, D ? =1.2x 10" 18 cm" 1 . Given the above pa- 
rameters, one would expect /^deflect ^438 AU, about where the 
southwestern lobe of the disk shows a warp of ^7°. 

We show a comparison between the F160W image of 
HD 32297 and our synthetic image in Figure [6] arbitrarily 
scaled to similar maximum brightnesses and using the same 
stretch. We also show a comparison between two contour 
maps of the surface brightness, using the same scaling. The 
model and the data agree qualitatively. 

4.1.2. An ISM Gas Interaction Model for HD 15115 

HD 15115 is an F2 star at 45 pc that has a disk with 
a highly asymmetric edge-on morphology; the Eastern lobe 
of the disk appears truncated relative to the Western lobe 
dKalas etal.ll2007t iDebes et al.ll2008al) . The proper motion 



vector of HD 15115 is southeast, with its eastern component 
almost twice as large as its southern component. In this case, 
we do not have a measure of the ISM cloud's radial veloc- 
ity, so we strictly consider the relative transverse motion in 
the sky between the cloud and HD 15115. We can convert 
this transverse motion to determine the relative velocity con- 
straints between the disk and the hypothesized ISM. The disk 
midplane is oriented with a PA of 279°. The transverse ve- 
locity vector of HD 15115 on the sky is at roughly the same 
PA as the disk. If we again assume a cloud with velocity 
of of 14 km/s on the sky, the minimum and maximum pos- 
sible relative transverse velocity is between 7 and 35 km/s 
if the cloud's velocity vector is parallel or anti-parallel to 
HD 15115's proper motion. The morphology of the disk is 
consistent with the midplane plowing headlong into an ISM 
flow. 

We assumed that HD 15115 had a stellar mass of 1.5 M© 
and that grains were launched from 35 AU, the rough in- 
ner radius of HD 151 15's disk as inferred from its IR e xce s s 
(IZuckerman & Songll2004t IWilliams & Andrewsll2006l) . Ex- 
perimenting with different models showed that D g = 1.3 x 10~ 18 
cm" 1 and a relative velocity vector with equal components 
along the midplane and along the observer's line of sight with 
a magnitude of 30 km/s best matched the observed data. 

Figure [7] shows a comparison of the F110W image of 
HD 15115 with our disk model on the same size scale. The 
eastern lobe (left side of Figure |7]i is contaminated by a 
strong residual that points roughly northeast in the Figure (see 
IDebes et al.l2008aT) . Otherwise, the figure shows that the mor- 
phology of the model matches the observation of HD 15115 
showing a significant truncation of one lobe of the disk. Ad- 
ditionally the figure shows that our models predict that the 
FWHM of the disk midplane should be wider on the truncated 
side. This aspect of the disk m orphology is not s ignificantly 
probed by current observations (Debes et al. 2008a). Another 
proposed scenario for this morphology is the close fly-by of 
HD 15115 by another star dKalas et al.ll2007l) . However, no 
definitive candidate s for a fly-by have currently been detected 
(IDebes et al.ll2008al) . 

4.1.3. An ISM Gas Interaction Model for HD 61005 

The disk observed around HD 61005, the "Moth", shows 
a peculiar disk morphology reminiscent of material being 
pushed up out of the disk midplane (H07). This shape might 
indicate disk interaction with the gas component of the ISM. 
The star itself is a late G8 dwarf at a distance of 34 pc, and 
H07 found that the surface brightness of the dust is consis- 
tent with dust grains roughly the size of the blowout radius, 
0.3 [im. H07 also determined from HD 61005's disk SED 
that the grains originate from radii ~10 AU from the star. 

The proper motion of HP 61005 is -56.09 m as in RA, 
and 74.53 mas in declination ( van Leeuwenll2007l) . If we ro- 
tate these velocities to a frame where the x-axis is parallel 
to the presumed major axis of the disk and convert to rel- 
ative velocities under the same assumptions we used above 
as for HD 15115, the relative velocities will range in magni- 
tude from 2 km/s to 30 km/s, with a vector primarily along 
the x-z plane at an inclination of 68° relative to the x-axis. 
If icy 0.3/im grains are the population we observe around 
HD 61005, we can derive D g under the assumption of an ISM 
gas density N(HI)=100 cm" 3 . These two assumptions result 
in D g = 3.3 x lO^cm" 1 

We experimented with various relative velocity vectors un- 
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til we found a qualitative match to the data. A velocity vec- 
tor that is nearly perpendicular to the disk midplane with a 
magnitude of 25 km/s matches the data best, along with an 
inclination between the observer and the disk of 20°. Figure 
[8] shows a comparison of the model and the data, at the same 
size scale and stretch. We also see good qualitative agreement 
with the morphology of this disk as well. 

5. DISCUSSION 

In this Section we discuss some details that our ISM inter- 
action model presently lacks. 

5.1. ISM Dust Grains 

For a disk interacting with an atomic ISM cloud with small 
grains, AC97 calculated an "avoidance radius", r av , where the 
radiation pressure of the central star was strong enough to 
scatter small grains away from a circumstellar disk. The dis- 
tance r av depends on j3, Voo, the initial velocity of a grain at 
large distances, and the stellar mass M+. 

2Q3-DGM* 
r m = 5 ■ (") 

Incoming ISM dust grains pile up at this radius and scatter 
away from the star, leaving a paraboloidal region with mini- 
mum radius r av that is free from sandblasting from small ISM 
grains. This boundary marks the separation between where 
grains from the ISM could have a significant effect and where 
the ISM dust does not impact the evolution of dust in a cir- 
cumstellar disk. 

The color asymmetries observed for HD 32297 could po- 
tentially be explained by this ISM grain sandblasting. With 
some assumptions about the ISM grains as well as the cloud 
that may be impacting HD 32297 's disk, we can estimate what 
r m is for HD 32297. AC97 calculated expected /3's for ISM 
dust for an A5 star such as HD 32297. They calculated a max- 
imum (3 for /3-Pictoris assuming silicate+carbonaceous grains 
of 19.2. HD 32297 has a similar luminosity to /3-Pictoris, so 
r av for HD 32297 would then be on order of between 133 AU 
and 1061 AU. The relative velocity of 15 km/s that we used 
for our gas interaction model of HD 32297 would correspond 
to an r av =245 AU. This is slightly larger than expected from 
the observations, and could be accounted for by a combina- 
tion of the gas interaction and sandblasting by ISM grains. 

AC97 calculated how ISM grain sandblasting would affect 
dust erosion rates in a disk. They found that there should 
be larger dust removal rates in the part of the disk directly 
impacting the ISM, which in turn causes a large asymmetry 
in dust production/erosion between regions that suffer direct 
sandblasting and those that avoid it. Additionally, there is an 
asymmetry between dust grains that are moving anti-parallel 
with the ISM flow and those that move away from the ISM 
flow, especially for grains in orbit around the star. Detailed 
modeling of the interaction with self-consistent grains of a 
particular composition and between the dusty ISM grains as 
well as the HI gas contained in the ISM cloud is left for a 
future paper (Debes, Kuchner, & Weinberger, 2009, in prep). 

5.2. Other Uncertainties in the ISM Interaction model 

While the interaction of dust with the surrounding ISM pro- 
vides an attractive explanation for the above disks, our models 
contain several uncertainties. For example, we have followed 
only a single size of dust grain in our interaction models, when 
clearly a size distribution of grains will be present in the disk. 



Furthermore, we have assumed a steady bulk flow to the ISM, 
even though one would expect turbulence on moderate size 
scales which could create further substructure. 

While we found specific values of D g , (3, and v re i that fit 
the observed morphology, any of these can be varied to bal- 
ance out a change in the other; we cannot derive unique fits 
to the data without further constraints on composition and 
v re i. Choosing a specific composition and grain type, such 
as spherical or fluffy particles would fix (3 and constrain D g . 
Modeling the morphology and surface brightness of a disk at 
multiple wavelengths may provide further constraint on a par- 
ticular gas interaction model. In some cases this effect should 
help to further constrain dust grain particles, since the exact 
morphology of an interaction depends on the dust properties. 

Physical processes, such as the presence of a stellar wind, 
can work to frustrate the interaction between the ISM and 
the circumstellar disk. For young solar-type stars, such as 
HD 61005, this may be important, while this may be less so 
for earlier type stars that have less magnetic activity. Regard- 
less, models of interactions between a dense ISM cloud and 
the sol ar wind can guide us to th e order of magnitude of this 
effect. iYeghikvan & Fahn d2003l) modeled the interaction be- 
tween an ISM cloud with N(HI)=100 cm" 3 and a solar-type 
wind. They found that when the sun was placed in such a 
dense cloud the heliopause shrunk to ~1 AU, well interior to 
the circumstellar region we observe for HD 61005. However, 
younger stars have stronger winds, which shield their debris 
disks out to hundreds of AU. 

Finally, there are large uncertainties in the velocities and 
densities of the modeled ISM clouds. Background stars near 
each of the above stars could provide a possible probe of the 
clouds and constrain their density/radial velocity by provid- 
ing lines of sight spatially coincident with the debris disks. 
Measurement of the transverse motions of the clouds would 
require a way of detecting the clouds directly at high spatial 
resolution, a difficult endeavor. 

6. CONCLUSIONS 

We have detected the HD 32297 debris disk at two addi- 
tional wavelengths, bringing to four the total number of wave- 
lengths in which this disk has been observed in scattered light. 
Both the morphology and the surface brightness of the disk as 
a function of wavelength hold useful clues to the nature of 
the dust in orbit around the star. The addition of more re- 
solved near-IR scattering images of the HD 32297 disk pro- 
vides measures of the morphology, surface brightness pro- 
files, and scattering efficiency of the dust at three separate 
wavelengths. 

We find that the morphology and asymmetries of the 
HD 32297 disk can qualitatively be explained by the inter- 
action of the disk's dust with a moderately dense ISM cloud. 
We developed a numerical model for the ineraction of a de- 
bris disk with ISM gas, and applied it to this system. We have 
also applied our ISM interaction model to two other disks, 
HD 15115 and HD 61005 and demonstrated that some warps 
and asymmetries in the outer regions of these disks can plau- 
sibly be explained by this mechanism. 

Future coronagraphic observations of HD 32297' debris 
disk, especially with STIS, will precisely measure the scatter- 
ing efficiency of the inner disk in visible wavelengths. More 
work needs to be done to understand precisely the effect of a 
debris disk's environment on how it appears. 

The authors would like to thank Dean Hines and Glenn 
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Table 1. Power-law indices to the Surface Brightness Profiles 



Filter Rbreak (AU) Inner Index (a) Outer Index (J3) 







Southwest Lobe 




F110W 


104±6 


-1.4±0.2 


-3.33±0.03 


F160W 


111±3 


-1.0±0.2 


-3.57±0.03 


F205W 


115±4 


-0.4±0.2 


-3.9±0.1 






Northeast Lobe 




F110W 


80±10 


-1±0.7 


-3.00±0.05 


F160W 


90±10 


-1.8±0.3 


-3.33±0.05 


F205W 


97±6 


-1.6±0.4 


-3.8±0.06 



Schneider for kindly providing an image of HD 61005 and 
useful discussions on both HD 61005 and HD 32297. We 
would also like to thank Carey Lisse for pointing out that stel- 
lar winds can affect dust/gas interactions and Chris Stark for 
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FIG. 1 . — Logarithmically scaled, final combined NICMOS coronagraphic images of the HD 32297 debris disk at Fl 10W (1.1/im, upper left), F160W (1.6/im, 
upper right), and F205W (2.05/jm, lower left). In the lower right corner of the figure is a logarithmically scaled color combined image of the disk, where blue is 
Fl 10W, green is F160W, and red is F205W. 
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FIG. 2. — Median measure of the position angle (PA) of each lobe of the HD 32297 debris disk as a function of distance from the star. A significant linear trend 
in the Southwest lobe is seen. 
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FIG. 3.— Surface brightness profiles of the Northeast and Southwest lobes of the HD 32297 debris disk in the F110W, F160W, and F205W filters. The 
F160W filters are scaled by a factor of 0.5, while the F205W filters are scaled by a factor of 0.125 for clarity. The surface brightnesses were measured in 3x3 
pixel apertures. A surface brightness asymmetry exists between roughly 1"5 to 3"in all three filters, with the asymmetry becoming more pronounced at longer 
wavelengths. 
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FIG. 4. — Scattering efficiencies at 3 different distances in the HD 32297 disk for the Northeastern lobe. 
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FIG. 5. — Same as Figure|4]but for the Southwestern lobe. 
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FIG. 7. — Comparison between an F110W image of HD 151 17's disk from Debes et al. 1 2008b) and a model of a disk interacting with a cloud of ISM gas with 
the parameters mentioned in Section l4. 1.21 
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FIG. 8. — Comparison between an Fl 10W i mage of HD 61005's disk (kindly provided by G. Schneider) and a model of a disk interacting with a cloud of ISM 
gas with the parameters mentioned in Section l4. 1.31 



